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In this paper we study thermalization in a strongly coupled system via AdS/CFT. 
Initially, the energy is injected into the system by turning on a spatially homogenous 
scalar source coupled to a marginal composite operator. The thermalization process 
is studied by numerically solving Einstein's equations coupled to a massless scalar 
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field in the Poincare patch of AdS§ . We define a thermalization time tx on the AdS 
side, which has an interpretation in terms of a spacelike Wilson loop (W(l ~ ^)) 
in CFT. Here T is the thermal equilibrium temperature. We study both cases with 
the source turned on in short(At < i) and long(Ai > i) durations. In the former 
case, the thermalization time = y < y and the coefficient gt — 0.73 in the limit 
At < ^P. In the latter case, we find double- and multiple-collapse solutions, which 
may be interpreted as the gravity duals of two- or multi-stage thermalization in CFT. 
In all the cases our results indicate that such a strongly coupled system thermalizes 
in a typical time scale tx — 
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I. INTRODUCTION 

Knowledge about thermalization in strongly coupled systems is important for us to un- 
derstand the results from relativistic heavy-ion collisions at RHIC and LHC and some non- 
equilibrium processes in the early universe. The AdS/CFT correspondence [l, 2] provides us 
with an elegant tool for studying the strongly coupled large Nc super Yang-Mills (CFT). 
The gravity dual of static plasma in 4-dimensional Minkowski space(M 4 ) is the AdS black 
brane(black hole) solution in the Poincare patch of AdS§. And studying gravitational col- 
lapse of matter fields in the bulk of AdS$ helps understand the thermalization process on 
the CFT side. 

To study thermalization by AdS / CFT, one needs first to choose some far-from-equilibrium 
states. Such states can be prepared by injecting energy into the CFT vacuum. This can be 
done by turning on sources for some boundary operators [3|. Examples for such boundary 
sources are the boundary metric coupled to the boundary stress energy tensor 4J or scalar 
sources coupled to marginal scalar composite operators [5j. Another kind of initial states are 
obtained by reconstructing some bulk metric based on the knowledge of the CFT datajsj]. 



One of this kind is the double shockwave metric |6Hl6l|. motivated by understanding ther- 
malization in heavy- ion collisions(see for a recent review). Of course one can also start 
from the gravity side by constructing consistent initial conditions for Einstein's equations. 



For this kind o 
cases and 119 



initial conditions, the interested reader is referred to 118] for boost-invariant 



23[ for the global patch of AdSd+i- Given initial conditions and suitable 
boundary conditions, thermalization in strongly coupled CFT may be studied by solving 
Einstein's equations in the bulk of AdS$. 



In this paper, we study thermalization of a spatially homogenous system in CFT 



24 



. The energy is injected into the system by turning on a spatially homogenous scalar 
source <fto{t) f° r a marginal scalar composite operator O. The cases of weak fields were 
first tackled by the authors of Ref. using perturbative techniques. We will deal with 
more general cases using a numerical method similar to that used in the global patch of 



AdSd+i 19H22J. In CFT the Lagrangian of the system in general takes the following form 

S = Soft + j d A x<j> (t)O{t) + j d 4 xA a J a , (1) 

thermalization probes 

where Soft is the CFT action, A$ and J a denote all the other external sources and cur- 
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rents(operators). On the gravity side, the bulk metric is determined by solving the grav- 
itational equations coupled to the massless scalar field 4> corresponding to cj) . And the 
back-reaction of all the other fields A a corresponding to the boundary sources Aq is as- 
sumed to be negligible. The expectation values of J a and their correlators can be calculated 
by solving the equations of motion of those weak fields in the bulk metric [21]. They are 
also essential for understanding the details of thermalization in a coordinate-independent 



way 



26]. 



In the introduction we summarize our results. The metric in Schwarzschild(Poincare) 
coordinates can be written in the form 

ds 2 = -1 (-fe~ 2S dt 2 + r l du 2 + dx 2 ) , (2) 
u 

where / and 5 are functions of t and u only. And the AdS black hole metric is given by 

u 4 

/wi = 1 j j and 5 = 0, (3) 

«o 

where uq = -\ and T is the thermal equilibrium(Hawking) temperature. The energy is 
injected into the system according toQ] 

T ( o 4) = (oWo, (4) 



where (Oj is the expectation value of O. The source 4>o(t) explicitly takes the following 
form 

Mt) = ~e- at \ (5) 

where e and a = are two parameters and At characterizes the duration of the source being 
turned on. However, our qualitative conclusions should hold regardless of the source's shape. 
The source induces an ingoing wave of the scalar field in the bulk, which will eventually 
collapse to form a black hole. We will not define the thermalization time tf in terms of the 



formation time of the apparent horizon because it is coordinate-dependent [29l. l30j . Instead, 
we define tr as the moment when the AdS black hole metric is established in the most part 
of the bulk that is causally connected to the boundary, can be interpreted as the scale- 



dependent thermalization time 26| defined by a spacelike Wilson loop (W{1)) with I ~ ^ in 



CFT. With such a definition, our results respect the following scaling invariance 

(a, e, tr, T) -»■ (A 2 a, A 2 e, t T /X, XT) (6) 
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FIG. 1. Narrow wave(At < ^). At t ~ At, the source on the boundary CFT is turned off and a 
narrow ingoing wave with a width An ~ At is induced near the boundary u = 0. The wave starts 
to propagate at the speed of light {du s /dt ~ 1) in the bulk and leaves behind a submanifold(n < u s ) 
equipped with the AdS black hole metric. It propagates more slowly in the deeper interior of AdS^. 
At t ~ du s /dt ~ and the AdS black hole metric is established in the most part of the bulk 
that is causally connected to the boundary. 



with A > 0, which gives 



9t 
T 



(7) 



g t is scale-dependent 



261 ] and in this paper we only discuss thermalization with / 



~ J. In 



the following, we will show that g t = 0(1) in this case. 

In the cases with ^ < 1, the boundary source induces a narrow wave with At = < ^ 
in the bulk. Our results are illustrated in Fig. [TJ We find that < 1 for all the narrow 
waves and g t ~ 0.73 if At < ^P. In CFT, the interpretation of our results is as follows: 
the system, after the source is turned off, thermalizes in a time scale tx ~ ^ in a top-down 



manner 
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QCD[31 



26] . This is in sharp contrast with bottom-up thermalization in perturbative 
32], in which soft gluons equilibrate more quickly than hard gluons. In CFT 
the interaction is equally efficient for soft and hard modes because of the vanishing beta 
function. At t ~ At, only time-like high momentum(u; ~ -h = y/a) modes are present in 
the system(see 33| for a similar interpretation in the case of a classical string and 34j for an 
R current). Those high u> modes split into low u modes very rapidly [34|. In the meanwhile 
the interaction is so efficient that the (remaining) high u modes equilibrate in a time ~ -. 
Such a splitting-equilibration continues from higher lo to lower oj modes. At t ~ ^ all the 
modes with uj > T achieve thermal equilibrium. 



In the cases with - > 1, the induced waves are broad and we do not evaluate the exact 
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FIG. 2. Two-stage and multi-stage thermalization. Fig. (a): The energy is injected into the CFT 
vacuum by the two pulses of 4>. If At > jr-, the system first thermalizes at a lower temperature 
Tl at t ~ 0. Then, the second pulse heats up the plasma to a higher temperature Tjj- Fig. (b) 
shows a source for a possible multi-stage thermalization. 



value of (7t in (EJ). The typical thermalization time ty shows up in a way illustrated in Fig. 
El The system is expected to thermalize in a time ty ~ -f- Therefore, if the source is turned 
on and off several times with period At ~ the system should achieve thermal equilibrium 
within each time interval At. We refer to such a thermalization pattern as a two-stage(Fig. 

or rnultiple-stage(Fig. EJ^b)) thermalization. This is the interpretation of the double- 
or multiple-collapse solutions that we find in the bulk. According to eq. (jlj), the energy is 
injected into the system by the two pulses of <f> as illustrated in Fig. EI a). After the source 
is turned off at t > 2At, the scalar field will eventually collapse. If - > 3.0, at t ~ the 
falling of the scalar field results in a submanifold defined by u < 0.98-Uo, which is equipped 
with the AdS black hole metric with the Hawking temperature T — T L . We interpret it as 
an intermediate thermal equilibrium state with T = Tl in CFT. And it is called the first 
collapse even though the trapped region has not formed yet (in Schwarzschild coordinates) at 
this time. Then the final collapse at t ~ 2At is naturally interpreted as the second stage of 
thermalization(heating-up process) in CFT. The criteria for obtaining such double-collapse 
solutions is At > If At > -r-, the system thermalizes in a time ~ At > ^_ simply 
because the thermal equilibrium is destroyed by the continuous injection of high u modes 
from top down before t = 0. A multiple-collapse solution can be defined in a similar way 
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and we find that the above criteria is also parametrically true for obtaining multiple-collapse 
solutions induced by a periodic source(see Fig. 12(b)). Therefore, we conclude that in such 
a strongly coupled system the typical thermalization time is tr ~ h- 

This paper is organized as follows. In Sec. [TTJ we derive the equations of motion for 
a massless scalar field coupled to gravity in the Poincare patch of AdS^. The numerical 
schemes in addition to initial conditions and boundary conditions are discussed in Sec. IIHI 
In this section, we also define the thermalization time £j>. Our numerical results are presented 
in Sec. |IVl In Sec. [V] we briefly conclude. In Appendix [A] we give the equations of motion 
in Eddington-Finkelstein coordinates. The details of our numerical methods are presented 
in Appendix [B] 



II. EINSTEIN-KLEIN-GORDON EQUATIONS 

On the AdS side, we need to calculate the back-reaction of a massless scalar field to 
the bulk geometry. The bulk action in AdSd+i corresponding to the first two terms on the 
right-hand side of is given by 



S = -4— ( / d d+l x^ {R-2K-2 («90) 2 } + 2 f d d x^K \ 

^ K d+1 U JdM J 



(8) 



where A = — -^jr- f° r AdSd+i, ^1 = AlT N 2 , L is a parameter of dimension of length, 7 is 
the induced metric on the boundary and K is the trace of the extrinsic curvature of the 
boundary. We need to solve the following Einstein-Klein-Gordon equations 

d a {V=99 ah d b <f>) = 0, (9) 

Rab — ^9abR ^T2 9ab = T a b, (10) 

where T a b is the stress tensor of the scalar field, which is given by 

T ab = 2d a ct>d b( t> - g ab (<90) 2 . (11) 



In the following, we take L = 1 and d = 4. 

In this paper, the scalar sources are assumed to be spatially homogeneous on the boundary 
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M 4 . In Schwarzschild coordinates, one needs to solve the following equations of motion 



V = u 



3 ( fe~*P 
u 3 



P = (fe- S V)' 



f = ^uf 2 e~ s VP, 

5' = \u (V 2 + P 2 ) , 

/ / = |«/(^ 2 + J P 2 ) + 1 (/-l) 



(12a) 

(12b) 
(12c) 

(12d) 

(12e) 



where the derivatives with respect to t and u are denoted respectively by overdots and 
primes, P = </)', V = / _1 eV and we take ([T2e]) constraint equation. If V and P vanish 
in a submanifold near the boundary of AdS§, the only solution to f!12dj) and (I12el) is 

,i 



/ = 1- 



u 



i • 



5 = 0, K = and P = 0, 



(13) 



which is Birkhoff 's theorem! 



in such a spatially homogeneous case in AdS 5 . The equations 



of motion above are invariant under the following scaling transformation 



x a ) ->■ (f) = (j)(\- l x a ), V ^ XV, P ->■ AP, 



x a ->■ x a 



Ax a , /(x a ) -)• /(x a ) = f(X x a ), 6(x a ) -> 5(x a ) = 5(A-V). (14) 



III. GRAVITATIONAL COLLAPSE OF MASSLESS SCALAR FIELDS 

In this paper, we study the response of the vacuum/plasma to an external scalar source 
in CFT. On the gravity side, the scalar source provides the boundary conditions and the 
vacuum AdS/the AdS black hole metric provides the initial conditions for solving eq. f|T2|) . 
We aim to see how a submanifold near the boundary of AdS^ equipped with the AdS black 
hole metric in ([3]) forms by gravitational collapse of a massless scalar field. 



A. Initial conditions 



Before the scalar source is turned on, the scalar field is assumed to vanish in the bulk. 
In this case, by Birkhoff 's theorem one has 

f bh = l-—, 5 = 0, V = and P = 0, (15) 
u 
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which is the AdS black hole metric in eq. fl3]), the gravity dual of static plasma. By 
holographic renormalization, the stress tensor of the boundary CFT is given by [35] 

T (4)^ = ^Ar c VT 4 diag{-l, 1/3, 1/3, 1/3}. (16) 

o 

Taking T — Y 0, one gets the vacuum AdS metric, dual to the CFT vacuum, as follows 

f vac = l, 5 = 0, V = and P = 0. (17) 
Both (TPTj) and (TT5l) will be used as initial conditions for solving ([12]) in this paper. 

B. Boundary conditions 

The boundary condition for solving (112dj) is given by 



5(t,0) = 0. (18) 

The boundary conditions for solving (I12al) and (112bl) are given by the scalar source in eq. 
§5§ in the boundary CFT, which is rewritten in the following form 

V(t, 0) = -2tee- at2 and P{t, oo) = 0. (19) 

The above boundary conditions give a unique solution in the bulk. Another solution can 
be obtained by replacing (a, e) with (A 2 a, A 2 e) in (fl9|) . These two solutions are related to 
each other by the scaling transformation in eq. f fT4|) . In the following, we denote such an 
equivalence briefly by 

(a,e) = (A 2 a,A 2 e). (20) 

As a result, we only need to study the dependence of solutions either on a or e(see Fig. [3] 
for an example). 

C. Numerical scheme 

The equations of motion in eq. f[T2"|) can be solved numerically either by the Chebyshev 
pseudo-spectral method jl5| or the finite difference method [361] (see Appendix [B] for a detailed 
description of our numerical methods). Given /, 5, V and P at t = t n , we first calculate 
V, P and / at the next time step t n+ i by solving f)12ap . (I12bj) and f)12cp . We use the third 
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FIG. 3. Usage of the scaling transformation in (|14p . Here, we show two numerical solutions with 
boundary conditions respectively given by (a, e) = (2.5,0.125) and (a, e) = (10,0.5). In these two 
figures the dashed curves show the metric functions / and 5 of the solution with (a, e) = (10,0.5) 
under the scaling transformation (A = 2). They are the same as those with (a, e) = (2.5,0.125). 
This example illustrates that for the equivalent solutions in (|20p we need only to calculate one of 
them and get the rest by the scaling transformation in (|14|) . 

order Adams-Bashforth method as our time-marching scheme. Then, 5 at t n+ i is obtained 
by solving eq. (112dj) . Given the initial conditions in ( ITT)) / ()15p . the bulk metric at late times 
can be calculated by repeating the above two steps. 

In numerical simulations, one needs only to study the evolution of the scalar field in the 
bulk region < u < u max with u max being some bulk cutoff. V = and P = in the bulk 
region in which the scalar field has not reached yet. From ( TT2")) . it is easy to show that / 
is independent of t and 5 is independent of u in this region. The geometry in this region 
can not influence the propagation of the scalar in the bulk. Therefore, one can arbitrarily 
choose u m ax to be any point in this region. In this case, the boundary condition at u = oo 
in ( TT9~)) is replaced by P(t,u max ) = 0. On the other hand, it helps to save computation 
time by choosing u max close to the deepest region that the scalar field can reach within 
the thermalization time t?. Using trial and error, we find that it is sufficient to choose 
u max ~ 2uo or u max ~ 3At respectively for narrow waves or broad waves discussed in the 
next section. 
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D. Thermalization time 

In this paper, we define a thermalization time tx by the first time when^l 

f{t T ,u) > f min and f(t T , ^min) fmini (21) 

ie black 



where f m in is chosen to be 0.01 for the narrow waves. A similar definition of t 
hole formation time is also used in the global patch by the authors of Ref.s 21 



22|. We 



also require that f(tT,u) at u < u m i n is that of the AdS black hole metric in f[To"|) . The 
goodness-of-fit is evaluated by 



E (Afc(«0//i - 1) 2 

, (22) 



n 

where fa is our numerical result at Ui and n is the number of the grid points Ui < u m i n . 
In the following we shall show that the energy density of the scalar field narrowly peaks 
around u = u min when u min > 0.98m . Before encountering the coordinate singularity at the 
apparent horizon(given by / = 0), we will refer to such a state as a collapse. 

The thermalization time has a coordinate-independent interpretation on the CFT 
side. It is the thermalization time defined by a non-local operator (0(1)) at scale I ~ ^. 
And the / m j n dependence of tr corresponds to the scale dependence of the thermalization 
time discussed in Ref. |26| . Here we only consider the rectangular spacelike Wilson loop 
(VK(Z))0, [^(see Ref. [26 1 for the discussion of other non-local operators). (W(l)) is a 
functional of spacelike geodesies which satisfy the following equations 

4 (^^) + (26fe- 2S d x t 2 - fe- 2S d x t 2 - f- 2 fd x u 2 ) = 0, (23) 

dx i!~^f) + ^ " 2dh ( 26 'f e ~ 25d ^ ~ f' e ~ 25d * t2 ~ r 2 f' d ^ 2 ) = °' ( 24 ) 

where G = a/1 — fe~ 2S d x t 2 + f~ 1 d x u 2 , d x = d/dx and the boundary conditions are given by 
u(— |) = = u(t) and t(— |) = t = t(|). As we shall see in the next section, at t ~ t T there 
is a sharp transition between the submanifold with vanishing V and P near the boundary 
and the rest of the bulk. In this submanifold, 5', 5 and / vanish. Therefore, it is a good 



1 As discussed below, tx is coordinate-independent but this simple form in the definition of tr is coordinate- 
dependent. See Appendix [A1 for the discussion in Eddington-Finkelstein coordinates. 



11 




FIG. 4. Submanifold probed by a spacelike Wilson loop. Fig. (a) shows u(0) as a function of I in the 
vacuum AdS(Vacuum) and the AdS black hole metric(AdS BH). Because of the scaling invariance 
in the vacuum AdS, u(0) is a linear function of I, i.e., u(0) ~ 0.835/(note that a logarithmic /-axis 
is used in Fig. (a)). Fig. (b) shows the spacelike geodesies u as a function of x with boundary 
conditions u(0) = 0.5,0.9 and 0.997. Here, we take no = = 1. 



approximation for us to solve the geodesic equation in the AdS black hole metric instead, 
which takes the following form 

d 2 2(l-2u* + u 8 + d x u 2 ) n 

-rr;U H ; tt = U, [25) 

dx 2 u(l-u A ) ' v ; 

where t(x) = t and we have taken uq = 1. Eq. (I25p can be easily solved numerically by 
the shooting method starting from u(x = 0) = u(0) and u'(0) = 0. As showed in Fig. 
S(a), u(0) is a monotonically increasing function of I. As a result, to probe high momentum 
modes(w ~ j 3> T) one only needs to know the bulk metric at u < I while for low momentum 
modes(u; ~ j ~ T) all the information of the bulk metric at u < u min ~ m is needed. In our 
definition of tr, u m i n = 0.997mo- Solving eq. ( |25l) with u(0) = u m i n gives Z = 3.16«o — ^r- 
Therefore, at t ~ the expectation values of (VF(7)) with Z < ^ all reduce to thermal 
results(see Fig. S(b)). In this paper, we only talk about thermalization in this sense and 



focus on the gravity side. A quantitative analysis like that in Ref. 26[ on the CFT side is 



also essential for obtaining more detailed information but we leave it to future studies. 
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FIG. 5. An ingoing narrow wave. This is a detailed example of the falling shell illustrated in Fig. 
[[J Fig. (a): f 2 (V 2 + P 2 ) and the metric function / at different times. Fig. (b): f 2 (V 2 + P 2 ) and 
^ = fe~ s , the speed of light-like geodesies, at different times. Here, (a, e) = (400,0.5). 

IV. RESULTS 

In this section, we discuss results with boundary sources of different values of (a, e) in 
( IT9|) . According to ( )20|) . we need only to study those parameter sets with different ratios of 
-. In the following, we first fix either e or a for convenience of our numerical calculation. 
Then, general results are obtained by the scaling transformation in (JUJ). We study both 
narrow {- < 1) and broad(^- > 1) waves. 

A. Narrow waves ( At = 4= < 1 < 1 

In the limit ^ < 1, the source induces a narrowly peaked ingoing wave propagating in 
the bulk of AdS^, which has been illustrated in Fig. [TJ Let us take for example a narrow 
wave with At = 0.05 and uq = 1.00. Fig. |5] shows the energy density of the scalar fielc0, / 
and the speed of light-like geodesies of the solution. At t ~ 2 At, the source on the boundary 
CFT is turned off and a narrow ingoing wave is induced near the boundary. The wave starts 
2 The energy density of the scalar field for an observer with n a = ^_ goo ffi) is P — n a n b T a i, — f (V 2 + P 2 ). 

In this paper, we use f 2 (V 2 + P 2 ) instead in order to show the energy density at different times in the 

same figure. 
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FIG. 6. The dependence of uq on a for e = 0.5(Fig. (a)) and / for a = 10 but different e(Fig. (b)). 
At a > 5.0, «o is almost independent of a. This fact and the scaling transformation in (j!4H allow 
us to conclude the relation in (|27p . which is verified in Fig. (b). 

to propagate with the group velocity du s /dt ~ 1 in the bulk. Its group velocity becomes 
slower in the deeper interior of AdS$. At t = 2.29 ~ ^p- the AdS black hole metric is 
established in the submanifold denned by u < 0.997mo and the speed of lightlike geodesies 
% = fe~ 5 = 0.0063 at u > Uq. This is interpreted as the gravity dual of the thermalized 
system defined by the spacelike Wilson loop (W(l — ^)) in CFT. At each time t > 2 At, 
the bulk metric is composed of three parts: (I) the AdS black hole metric behind the wave, 
(II) the Vacuum AdS metric with a time dilation, i.e., 5 > 0, before the wave and (III) the 
transition over a region An < AAt between (I) and (II) across the wave. Therefore, in this 
case the system in the boundary CFT thermalizes in a top-down manner [26 1. Fig. E(b) 
shows that the speed of light becomes slower at later times in region (II). 

Let us first fix e = 0.5. As showed in Fig. E^a), uq only weakly depends on a as a > 5. 
Using the scaling transformation in eq. (|14p and neglecting such a weak dependence on a, 
one can conclude that as long as the condition 

- < 0.1 (26) 

(X 

holds, Uq is (almost) independent of a and proportional to e~^, that is, 



u = -2= ~ 0.699e^5 ; OI) T = ~ 0.455e5, (27) 

v2e Tiun ' 
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FIG. 7. tx and uq for e = 0.5(Fig. (a)) and g t as a function of ^-(Fig. (b)). 



where Mq ' 5 '* is our numerical result with (a, e) 



i( a5 ),0.5) and a^ ^ denotes the value of 
a with fixe e = 0.5. The above relation T oc is also obtained by perturbative techniques 
in the limit - <C 1 in Ref. eJ. Therefore, the perturbative techniques should be applicable 
when eq. ( |26i) is satisfied. The condition ( 1261) can be further rewritten in the following form 

(28) 



1 . 0.14 

At EE — < — - . 

\Ja T 



As a verification of ([27)1 . Fig. Mjo) shows the metric function / with e = 0.1,0.05 and 0.01 
and fixed a = 10 at two different times. In region (I) they all agree with the AdS black hole 
metric with uq given by ( 1271) . 

Fig. UJa) shows the thermalization time tx for e = 0.5 and different values of a. Again, 
the scaling transformation in f JT4j) allows us to make conclusions for the cases with (a, e) = 
(AV ' 5 \0.5A 2 ) S (a( a5 \0.5), that is, 



,(0.5) 



7TU, 



(0.5) T 



T 



(29) 



where g t ~ 0.73 for ^ > 200 or, equivalently, At < ^ and g t < 1.0 for - < 0.5(see 
Fig. UJb)). Note that the thermalization time tx is physically different from the coordinate- 
dependent formation time of the black brane tsH ~ At in ingoing Eddington-Finkelstein 
coordinates reported in Ref. 5[ . We discuss this in detail in Appendix |A] 

In the limit a — > oo, one may make a simple estimate as follows: the source induces a 
three-peak wave in the bulk as showed in Fig. [5j The last (leftmost) peak does not contribute 
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so much to the bending of the spacetime but does determine the location of the turning point 
between (I) and (II). Since the leftmost piece of this peak propagates in an approximate AdS 
black hole background, the thermalization time tx may be estimated by 

nin du 1 / U m i n , Un + U m i n \ 0.64 . . 

— = -u 2 arctan — + log — ~ 2.02u ~ — . 30 

Jbh 4 \ M U -U min J T 

At worst, this is the lower bound for tx, which guarantees that the system thermalizes 
without violating causality. We take this as another justification for defining tx by eq. ( 12~TT) . 
Taking a = 400, our numerical calculation gives tx — 2.29 and the estimated thermalization 
time is given by (At T + 2 At) ~ 2.12. 

B. Broad waves(Ai > ±): ^ > 1 

In this subsection we discuss solutions with - > 1, or, At > ^ by eq. ( 12"B"j) . As illustrated 
in Fig. EJ^a), the energy is injected into the CFT vacuum by the two pulses of 0. Based on 
the discussion in the previous subsection, the system is expected to achieve thermalization 
in a time scale tx ~ ^ • Therefore, if At > ^, one can expect a two-stage thermalization. 
Similarly, if <p is periodic with period At > ^, a multi-stage thermalization, as illustrated 
in Fig. Mb), may also be expected. We shall show how the results on the gravity side live 
up to such an expectation. 

In the cases with — > 1, it is more convenient to study solutions with a fixed at a smaller 
value. We choose a = 0.1(At = 3.16). Solutions with 1.0 > e > 0.05 are obtained by 
numerically solving ffl2|) . We observe a transition from a single collapse to a double collapse 
in those solutions at e ~ (0.2 — 0.3)[]. Fig. |8] shows the metric function / of a single-collapse 
solution(e = 0.1) and a double-collapse solution(e = 0.3). In the single-collapse solution the 
scalar field collapses once at u$ = 1.76 and t = 6.61 ~ 2At. In contrast, in the double- 
collapse solution the scalar field collapses twice respectively at u ~ ul = 0.73 and t ~ Z 
and at u ~ u H = 0.48 and t = 5.47 ~ 2At. For the first collapse, ^ = nu L = 2.29 < At 

3 Here, the first collapse at t = is determined by the following requirements: (a) f(0,u) > f m in for all u, 
(b) f(0,u m i n ) = fmin and (c) / is that of the AdS black hole metric at u < u m i n . The exact transition 
value of e depends on the choice of f m in or, equivalently, the scale I for (W(l)). For the solution with 
- = 3.0 we have / ~ or, equivalently, u m i n ~ 0.98uo f° r the first collapse. At t = 0.61, min{/} ~ 0.01. 
However, / is not that of the AdS black hole metric near the boundary. Therefore, we do not interpret it 

as being in thermal equilibrium. 

4 Here, ml is obtained by the least-squared fit of fbh with uq — ul to our numerical results for u < u m i n - 
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(a) (b) 



FIG. 8. Single-collapse(Fig. (a)) and double-collapse(Fig. (b)) solutions. In both figures we show 
/ as a function of u from t = — 3.0(top right) to 7.0(bottom left) with time steps equal to 1. In 
Fig. (a), / at t = is far from that of the AdS black hole metric. In Fig. (b), / at t = is that of 
the AdS black hole metric at u < u m i n = 0.98uo- 

while for the single-collapse solution 2 At > ^ = 5.53 > At. The scalar fields induced by 
the source with e > 0.3 all undergo such a double collapse. 

Let us understand better the double-collapse solution. In Fig. Owe plot the energy density 
of the scalar field, / and the speed of lightlike geodesics(du / dt = fe~ 5 ) of the solution with 
e = 0.3. At t > 5.47, two peaks are observed in the energy density. They correspond to 
the two collapses at u ~ ul and u ~ uh- Also the speed of light drops dramatically in the 
bulk region u > uh- The peak at u ~ Ul is resulted from the first collapse at t ~ when 
the speed of light starts to drop below 10~ 6 at u > u^. From t — to 5.47, the scalar field 
mainly accretes at u ~ %, which only significantly modifies / at u < ul and results in the 
second peak of the energy density. 

As another support for regarding the state at t = as being in (approximate) thermal 
equilibrium, we study the response of static plasma with T = T L to the scalar source 
which is turned on only after t = 0. On the AdS side, we solve eq. f lT2"|) using the AdS 
black hole metric in (TT5"j) as the initial conditions. Here, we only discuss the solution with 
T = Tl = 0.73 and (a, e) = (0.1,0.3). Fig. [TOl shows the metric function / at different 
times and the energy density at t = 6.0 of this solution. Their small difference from those 
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FIG. 9. The energy density of the scalar fleld(/ 2 (y 2 + P 2 )), / and the speed of light {du/dt = fe s ) 
for the double-collapse solution with (a,e) = (0.1,0.3). 

of the double-collapse solution with (a, e) = (0.1,0.3) justifies that an intermediate thermal 
equilibration is established at t ~ 0. 

Boundary conditions with - > 3 all give such double-collapse solution^ by the scaling 
transformation in f|T4l) . In general we conclude that the criteria for the system to achieve 
thermalization induced by the first pulse of 4> at t ~ is At > jr. This is exactly what we 
expected on the CFT side. The intermediate thermalization temperature Tl for arbitrary 
(a, e) is given by 

... 1 v/Ifti 



where (j^-) , showed in Fig. HU is our numerical results with a = 0.1. At larger e " ! 
scales according to the following power law 

4 ai) (e) -O.lOOe^ (32) 
5 We have not calculated solutions with even larger - > 10 because it is numerically too time-consuming. 
However, based on the physical argument given in the paper, it is reasonable to expect that the conclusion 
should also apply to the cases with | > 10. 



> L 
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Double Collapse ■ 
Black Hole 



0.2 



0.4 



0.6 



0.8 



FIG. 10. Collapse of a massless scalar field on the AdS black hole background. In both figures, 
the dotted curves (Black Hole) show the results of turning on the scalar source at t > on the AdS 
black hole background. The solid curves(Double Collapse) show the results of the double-collapse 
solution with (a, e) = (0.1,0.3). In the left figure, / is showed as a function of u from t = 0(top 
right) to 6.0(bottom left) with time steps equal to 1. 




o.oi 



FIG. 11. function of e/10a. 



with j3 ~ 1.74. The second collapse(the heating-up process) can also be studied by turning 
on the source at t > and starting from the initial conditions in ( IT5|) with uq = ul = -^r-, 
as we did for the case with (a, e) = (0.1,0.3). 

The qualitative conclusion above is independent of the source's shape. As a justification, 
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t=ir, u =1.52 
t=2.0ir, u =1.30 
t=3.0TT, u =1.17 
t=4.07r, u =1.06 
t=5.0 7t, u =0.99 
t=6.07T, u =0.92 
t=7.07T, u =0.87 
t=8.0rr, u =0.82 
t=9.07T, u =0.78 
t = 10.0 71, u„=0.74 



FIG. 12. A multiple-collapse solution. Fig. (a) shows / and f^h at different times from t = 7r(right) 
to t = lOvT(left) with time intervals equal to tt. Here uq in is obtained by performing least- 
square fit of the AdS black hole metric to our results at u < 0.505 with a ~ 10 -4 . Fig. (b) shows 
the energy density of the scalar field (/ (V 2 + P 2 )) at t = Wtt. Each spike here is the result of the 
collapse around t = it, 2tt, ■ ■ -, Wtt. 

Fig. [T2] shows a multiple-collapse solution. The periodic source is given by 

0o = esin(t)0(t) (33) 

with e = 0.3. Here, one naturally takes At = tt. Fig. [T2T a) shows / at t = it, 2tt, ■ ■ •, 10n. 
At t = 7r, / is differento from that of the AdS black hole metric at u > u m i n = 0.93mo 
and i = ttuq = 1.52/T > At. In contrast, at t = IOvt, / is that of the AdS black hole 
metric at u < 0.99mo an d -f = 0.74/T. The energy density of the scalar field at t — 10tt 
is showed in Fig. IT2T b). Each spike in the figure is the result of the collapse around 
t = tt, 2vT, • • •, 10/T. Therefore, the parametric criteria for a system driven by a periodic 
source to achieve thermal equilibrium within each time interval At is also At > ^ with T 
being the intermediate thermal equilibrium temperature. 

V. DISCUSSIONS 



In this paper, we focus on the bulk metric resulted from gravitational collapse of a 
massless scalar field in the Poincare patch of AdS^. Mostly, we aim to understand the 



6 Here, we define w mjrl by 



fbh(Ur, 



- 1 



= 0.1 with f m in being our numerical result at u = u r 
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typical thermalization time scale of its CFT dual. We find that thermalization in such a 
strongly coupled system is rapid in the sense that 

* - ^, (34) 

where the coefficient of 0(1) depends on the scale I ~ A of nonlocal operators and on 
the boundary source's shape. Such a rapid thermalization time seems to be typical of the 
strongly coupled CFT^, 38 1. 



We leave many unanswered but intriguing questions for future studies. There are still 
more details about such a thermalization process that can be understood only by evaluating 



non-local operators 



33 
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40| . Moreover, in this paper we only 



26l | or other relevant probes 
study gravitational collapse of massless scalars. It is interesting to know what is the typical 
thermalization time tx in the cases of massive and tachyonic scalars, which are respectively 
dual to irrelevant and relevant operators in the boundary CFtQ. 

ACKNOWLEDGEMENTS 

The author would like to thank R. Baier, Y. Kovchegov, P. Romatschke, S. Strieker and 
A. Vuorinen for reading this manuscript and providing illuminating comments. This work 
is supported by the Humboldt foundation through its Sofja Kovalevskaja program. 

Appendix A: Equations of motion in Eddington-Finkelstein coordinates 

In outgoing Eddington-Finkelstein coordinates the metric can be written in the following 
from 

ds 2 = — (-fe~ 2S dv 2 + 2e~ s dvdu + dx 2 ) , (Al) 
u v ' 

where / and 5 are functions of v and u only. We only need to solve the following equations 
of motion 

f = ^uV(fP + e 5 V), P = V, 



5' = -uP 2 , V = --- -v* -fe~ d P , (A2) 
3 2 u 2 \u 6 J 

where the derivatives with respect to v and u are denoted respectively by overdots and 
primes, 

P = (f/ and V = <p. (A3) 
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II: Vacuum AdS 
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FIG. 13. Coordinate dependence of the boundary time. Fig. (a): An event occurs at u = u s and 
at t in Schwarzschild coordinates. Atf, is the time difference from t on the boundary in outgoing 
Eddington-Finkelstein coordinates. Fig. (b): Results of / calculated in these two coordinate 
systems. Here (a, e) = (10, 1), u s = u m i n = 0.685 and uq = 0.697. The event in this example is 
that / drops to f m i n = 0.07 at u m i n for the first time. 



One constraint equation is given by 



f = \uP*f + -U-l). 

6 U 



(A4) 



The simple form of the definition of tx in eq. ( 121]) is due to our choice of coordinates. If 
we naively define tx by the same definition in ( 1211) in both coordinate systems, the thermal- 
ization time £j> calculated in Eddington-Finkelstein coordinates in (|A1[) is larger than that 
in Schwarzschild coordinates in ()2]) with a time difference given by 

At fe = / — = -uo 2 arctan h log , (A5) 

JO J 4 V U Uo-UminJ 

where we have taken f — 1 — m 4 /mq. The reason for this is illustrated in Fig. flUT a). Fig. 
[151(b) shows one example with u min = 0.685 and u = 0.697. In this case, eq. ( 1A5j) givesll 



At b ~ 1.10, 



(A6) 



which is exactly what we find in our numerical results. Therefore, to get the same tr in 
Eddington-Finkelstein coordinates one has to solve the spacelike geodesic equation (76) in 



7 Here, we choose f(t or v,u m i n ) — f m %n with f m i n — 0.07 for saving computation time. 



22 



Ref. 26j, which complicates our numerical calculations. Similarly, in ingoing Eddington- 
Finkelstein coordinates |^| such as that used in Ref. jsj one can expect a shorter formation 
time of the black brane. 



Appendix B: Numerical methods 



Let us write dt and du for the line spacings and denote the approximate value of any 
quantity h(t,u) at the grid point (t, u) = {ti,uj) by 



h l - ~ 



h(ti,Uj) 



(Bl) 



where t t =t Q + idt with i — 0, 1, . . ., to is the initial time and Uj is given below. We use the 
following third-order Adams-Bashforth method as our time marching scheme 



h\ H 

J 12 



23h) - 16^ 



i-l 



5h 



i-2 



(B2) 



In the following we discuss the choices of the grid points Uj in both the Chebyshev pseudo 
spectral method and the finite difference method. We will only focus on solving eq. ([12] 
but eq. (1A2I) can be solved in a similar way. 



1. The Chebyshev pseudo-spectral method 



In this case, the grid points Uj are chosen to be 

— J and j = 0, l,...,n. (B3) 

Accordingly, the differentiation matrix 

D n+ i = 2d n+1 /u max , (B4) 

where the elements of the (n + 1)— by— (n + 1) matrix d n+ i are given by( Ijj] and references 
therein) 



(B5) 



2ajti i=j = 

< i = j < n 

1 — j — n 

However, in ingoing Eddington-Finkelstein coordinates, say, ds = -^ (—fe~ 2S dv 2 — 2e~ s dvdu + dx 2 ) it is 
difficult to provide initial conditions for numerical simulations in our case. 



2n 2 + l 
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with Co = c n = 2 and Cj = 1 otherwise. The difference equations are obtained by replacing 
(. . .)' with D n+1 (. . .) in f|T2|) . Given V, P, f and 5 at tj and the previous time steps, all 
VJ +1 and -Pj +1 are calculated by solving the discretized versions of (pa) and ffT2bl) except 
V^ +1 and P i+1 , which are given by the boundary condition in ( TT9l) . i.e., 

V* +1 = -2t i+l ee-<+^ and P m = 0. (B6) 

2. The finite difference method 

In this case, the grid points Uj are chosen to be equally spaced, that is, 

Uj = jdu + u min with j — 0,1, ... ,n and du = [u max - u min ) / n. (B7) 
We take u m i n = 0.005, which is introduced to avoid the numerically singular behavior near 



u = 22J. The derivatives with respect to u in eq.s fl 1 2aj) to f!12dl) are replaced by finite 
differences 

Vj = —fie-Qpj + (.// "P - /; ;< ' P ; ) /du, (B8a) 

Pj = (/;.,' v; +1 - ./> 17) /d«, (B8b) 
/; ^V']- ( B8c ) 

^ " 5U ~ |«i (V7 2 + if) • (B8d) 



rfw 3 

With the initial conditions in ( TT7|) / (1T5]) and the boundary conditions 

V^ +1 = -2t m ee- a *?+i, P^ +1 = 0, and 5 l +1 = 0, (B9) 
one can calculate the late-time geometry by (JB8]) and (1B2[) . 
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